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The classical Heisenberg antiferromagnet with uniaxial exchange anisotropy, the XXZ model, in a 
magnetic field on a simple cubic lattice is studied with the help of extensive Monte Carlo simulations. 
Analyzing, especially, various staggered susceptibilities and Binder cumulants, we present clear 
evidence for the meeting point of the antiferromagnetic, spin-flop, and paramagnetic phases being 
a bicritical point with Heisenberg symmetry. Results are compared to previous predictions based 
on various theoretical approaches. 
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Uniaxially anisotropic Heisenberg antiferromagnets in 
a magnetic field have been studied quite extensively in 
the past, both experimentally and theoretically [J 
Usually, they display, at low temperatures and fields, 
the antiferromagnetic phase and, when increasing the 
field, the spin-flop phase. A prototypical model de- 
scribing these phases as well as, possibly, multicritical 
points, is the Heisenberg model with a uniaxial exchange 
anisotropy, the XXZ model 

n^jJ2[^is:s^ + sfs]) + s^s;] ~ hJ^s! (i) 

where J(> 0) is the exchange coupling between classi- 
cal spins, , S'^^.-^ , S'f(^.) ) , of length one at neighboring 
sites, i and j, on a simple cubic lattice, A is the uniaxial 
exchange anisotropy, 1 > A > 0, and H is the applied 
magnetic field along the easy axis, the z~axis. The phase 
diagram of the model has been investigated already sev- 
eral years ago, using, among others, mean-field theory 
3 , Monte Carlo (MC) simulations [3] , and high temper- 
ature series expansions The transition between the 
antiferromagnetic (AF) and spin-flop (SF) phases seems 
to be of first order, while the boundaries of the paramag- 
netic (P) phase to the AF and SF phases are believed to 
be continuous transitions in the Ising and XY universal- 
ity classes, respectively. Moreover, based on renormaliza- 
tion group analyses in one- loop-order, a bicritical point 
in the Heisenberg universality class had been proposed, 
at which the three different phases meet Q . This sce- 
nario has been questioned on the basis of renormalization 
group calculations in high-loop-order Q, where the bi- 
critical point has been argued to be unstable against a 
tetracritical point 9], which, in turn, may be unstable 
towards transitions of first order in the vicinity of the 
meeting point of the three phases. However, a subse- 
quent renormalization group analysis in two-loop-order 
[10] suggests that a bicritical point in the Heisenberg uni- 
versality class can not be excluded. 

As has been noted quite recently not only AF 
and SF phases, but also biconical (BC) structures [7] 
may play an important role in the XXZ model. Indeed, 
such BC structures are degenerate ground states at the 
critical field separating AF and SF configurations at zero 




^■^1 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.1 
kJ/J 

FIG. 1. Phase diagram of the XXZ model on a cubic lattice 
with exchange anisotropy A = 0.8 in the vicinity of the AF- 
SF-P point, as obtained from the present MC simulations. 



temperature. For the XXZ model on a square lattice, 
these degenerate BC fluctuations lead to a narrow disor- 
dered phase intervening between the AF and SF phases 
at low temperatures, giving, presumably, rise to a 'hid- 
den tetracritical point' |ll[Jl3| at zero temperature. A 
recent Monte Carlo study 13] for the XXZ antiferromag- 
net on a simple cubic lattice showed that biconical struc- 
tures, arising from the degenerate ground states, also oc- 
cur at low temperatures close to the transition between 
the AF and SF phases. But they do not destroy the di- 
rect transition of first order between these two phases, 
in accordance with the behavior predicted by mean-field 
and other, more reliable theories. 

In that recent MC study [11] of the three-dimensional 
XXZ model, Eq. (1), at fixed anisotropy, A — 0.8, con- 
tinuous phase boundaries of Ising type, for the AF-P 
transition, and of XY type, SF-P, as well as the transi- 
tion line of first order, AF-SF, have been identified. The 
position of the meeting (or AF-SF-P) point has been es- 
timated, keT/J = 1.025 ± 0.015 and H/J = 3.9 ± 0.03, 
without determining critical properties of the AF-SF-P 
point. The aim of the present paper is to deal with this 
intriguing aspect, using extensive MC simulations. 

Here, the standard Metropolis algorithm fl3] with sin- 
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FIG. 2. Isotropic Binder cumulant V^^^ at ksT / J = 1.0225, 
varying the field, for lattice sizes L = 16 (circles), 24 
(squares), and 32 (diamonds). 



FIG. 3. As for Fig. 2, at ksT/J = 1.025, including {V'^'y 
as dashed line [H, [l^ . 



gle spin-flips is applied. Employing full periodic bound- 
ary conditions, lattices with sites, are considered. L 
ranges from 4 to 40, with the main focus on the sizes 
L= 8, 16, 24, and 32, to study systematically finite-size 
effects. MC runs with, at least, 10^ MC steps per site 
for the larger lattices, are performed. Error bars are es- 
timated by averaging over, at least, three independent 
realizations. In this way, data of the desired accuracy 
are obtained, and there is no need to use other, perhaps, 
more powerful MC algorithms. 

To map the phase diagram, we study, especially, quan- 
tities related to the Ising, XY, and Heisenberg order pa- 
rameters. In particular, we monitor, apart from the abso- 
lute values of the staggered longitudinal, mf^, transverse, 
m^f, and isotropic, m^f^, magnetizations, the corre- 
sponding staggered susceptibilities, Xst^ Y»f i and Xsf as 
well as the corresponding Binder cumulants [la] , t/^ , U^^ , 
and U^^^ . To our knowledge, the isotropic quantities, be- 
ing crucial to identify a bicritical point with Heisenberg 
symmetry, had not been included in previous MC simu- 
lations. 

The phase diagram, in the {ksT / J, H/ J)-plane, in the 
vicinity of the AF-SF-P point is depicted in Fig.l. As 
before, we set A= 0.8. The shown estimates for the tran- 
sition lines are based on usual finite-size analyses |16|] . In 
particular, to determine the AF-P boundary, with the 
transition being in the Ising universality class, the size 
dependence of the height and position of the maximum 
in the longitudinal staggered susceptibility, xlt^ allows 
for reliable estimates. U^^ turns out to be very useful in 
estimating the SF-P transition line. There one observes 
rather small finite-size correction terms to the critical 
Binder cumulant for cubic lattices in the XY universality 
class [13, (C/^f)*= 0.586.., with the critical Binder cu- 
mulant being the cumulant at the transition in the ther- 
modynamic limit, L — cx). Finally, the AF-SF phase 
boundary of first order is readily identified from the lo- 
cation of the maxima in x^st and x^st ■ 

The resulting phase diagram confirms and refines pre- 
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FIG. 4. As for Fig. 3, at fcsT/J = 1.03. 



vious, independent MC findings [Ij]- Note that the 
present results allow to locate the AF-SF-P point accu- 
rately, keT/J = 1.025 ± 0.0025 and H/J = 3.89 ± 0.01, 
reducing appreciably the error bars obtained in the pre- 
vious MC study, as mentioned above. 

Most importantly, critical properties of the AF-SF- 
P point are studied. In case of a bicritical point 
with Heisenberg symmetry, the critical Binder cumulant 
^jjxyzy expected [H, 111 to acquire the value 0.620 
(±0.003), using periodic boundary conditions for lattices 
of cubic shape. Note that, in general, the critical Binder 
cumulant, in a given universality class, may depend on 
boundary condition, system shape, and lattice anisotropy 
of the interactions 15, 20, 21]. 

In Figs. 2, 3, and 4 the size dependent isotropic Binder 
cumulant U^^^' is displayed near the AF-SF-P point, 
at three fixed temperatures, ksT/J— 1.0225, 1.025, and 
1.03, varying the field to cross the boundary of the AF 
phase, see also Fig. 1. We may distinguish three different 
scenarios. At the lowest temperature, ksT/J = 1.0225, 
see Fig. 2, the cumulant increases with larger system 
sizes at all fields, so that there are no intersection points 
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FIG. 5. Staggered isotropic susceptibility Xsf^ ksT / J — 
1.025 for L= 16 (circles), 24 (squares), and 32 (diamonds). 
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FIG. 6. Effective critical exponent {'y /u)''^^ (Le) of the stag- 
gered isotropic susceptibility, see Eq. (2), at ksT/J = 1.025 
and H/J= 3.865 (circles), 3.875 (squares), 3.8875 (diamonds), 
and 3.90 (triangles). The dashed line marks the asymptotic 
critical exponent in the Heisenberg universality class • 



for cumulants of different sizes. This behavior is in ac- 
cordance with being at a temperature below that of the 
AF-SF-P point. Indeed, W^^^ tends to 2/3 in the AF 
and SF phases in the thermodynamic limit, and there 
is no indication of a transition of Heisenberg symme- 
try. At ksT/J = 1.025, see Fig. 3, all intersection 
points of the cumulants, for lattices of sizes L= 16, 24, 
and 32, occur closely to the critical Heisenberg value, 
^jjxyzy 0.620. This fact may be interpreted as ev- 
idence for being in the immediate vicinity of a bicriti- 
cal point of Heisenberg symmetry. Actually, moving to 



higher temperatures, another scenario shows up. There 
are still intersection points of the cumulants for differ- 
ent lattice sizes, but they shift for larger lattices to lower 
and lower values below that of the critical cumulant in 
the Heisenberg universality class. This trend is already 
seen for ksT/ J — 1.0275, and it is quite pronounced at 
kBT/J = 1.03, as depicted in Fig. 4. In any event, the 
observations on the isotropic Binder cumulant are com- 
pletely consistent with the existence of a bicritical point 
with Heisenberg symmetry at fc^T/J = 1.025 ± 0.0025. 

This suggestion is supported by the behavior of the 
isotropic staggered susceptibility Xst^ ■ Simulation data 
at fixed ksT/J = 1.025 are shown in Fig. 5. At first 
sight, there seems to be no hint for criticality of Heisen- 
berg type, because Xs^^ as a function of field exhibits no 
maximum. However, additional information is provided 
by the standard discrete effective critical exponent, at 
given temperature and field, 

{l/i^r^iL,) = HxT{Li)/xT{L2))/HLi/L2) (2) 

for successive sizes Li > L2 and with the effective length 
Lg — \/LiL2. In case of critical behavior in the Heisen- 
berg universality class, the effective exponent would ap- 
proach, for Le — >■ 00, the known (2^ asymptotic value 
i'^/vY^^-— 1.96... In fact, as shown in Fig. 6, the ef- 
fective exponent seems to approach closely that value 
at H/J = 3.8875 ± 0.0125. Obviously, at those fields 
one is in a (crossover) critical region governed by the 
Heisenberg fixed point, in accordance with the findings 
for the isotropic Binder cumulant. The suggestion is 
corroborated by the behavior of the staggered longitudi- 
nal and transverse susceptibilities. Both quantities show 
pronounced maxima, when varying the field at k^T j J— 
1.025. The size dependence of the peak position al- 
lows us to estimate the critical field. Within the error 
bars, both susceptibilities lead to the same critical field, 
Hj J = 3.895 ± 0.005, indicating, again, the closeness of 
the bicritical point. 

It should be mentioned that we did, in addition, a 
preliminary MC study at a rather strong anisotropy, 
A — 0.2. The AF-SF-P point is estimated to be at about 
ksTlJ = 0.2095 ± 0.0005 and HjJ = 5.825 ± 0.0025. In 
the immediate vicinity of that point, the isotropic Binder 
cumulant shows similar features to the ones discussed 
above. We tend to believe that the AF-SF-P point in 
the XXZ Heisenberg antiferromagnet on a simple cubic 
lattice is generically a bicritical point with Heisenberg 
symmetry. 
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